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We  have  completed  computer  simulations  of  the  inner  magnetosphere,  for  the 
substorm-type  event  of  19  September  1976.  Several  computer  runs  have  been  per¬ 
formed,  with  various  sets  of  input  parameters.  Data  from  the  Air  Force  S3-2 
satellite  and  ground  magnetometers  were  used  for  model  inputs.  Model  predictions 
vere  analyzed  and  exhaustively  compared  with  data  from  the  S3-2  satellite,  the 
AFGL  magnetometer  chain  and  many  other  sources.  Agreement  between  theory  and 
data  was  very  satisfactory  in  most  respects.  The  model  injected  a  realistic  * 

;ring  cu’-rent  during  the  substorm.  Many  features  of  various  type_s_of  observations ,jv 


^  /I'SIFICA'ION  .  i  n,B  _ _ _ _ _ _  — - - 

*were  "correctly  predTcTe3“by  The  model ,  which  was  then  used  to  suggest  physical  in¬ 
terpretations  of  the  features.  The  predicted  latitudinal  distribution  of  region-  ! 

2  Birkeland  currents  disagrees  with  S3-2  observations  in  a  way  that  appears  to  be  , 
independent  of  details  of  the  model,  input  inaccuracies,  etc.  Some  implications  j 
of  this  discrepancy  for  the  theory  of  the  plasma  sheet  and  magnetospheric  sub¬ 
storms  have  been  worked  out.  The  computed  total  Joule  heating  rate  in  our  I 

modeling  region  is  larger  than  would  be  estimated  by  a  standard  method,  and  the 
time-integrated  Joule  heating  rate  is  about  twice  the  change  in  model  ring-current 
energy. 

Recently,  our  program  has  been  converted  for  use  in  modeling  a  magnetic  storm 
and  one  run  has  been  made  through  the  first  few  hours  of  the  magnetic  storm  of  29  i 
July  1977.  A  realistic  model  ring  current  was  injected  in  the  first  few  hours  i 
after  the  sudden  commencement.  <e - 

A  large  computer  program  ha's  been  developed  that  integrates  the  Biot-Savart 
law  over  the  current  system  predicted  by  our  overall  event  simulation,  to  produce 
"theoretical  magnetograms"  for  the  event.  These  theoretical  magnetograms  have 
been  computed  for  the  AFGL  northern-U.S.  chain,  as  well  as  various  high-and  low- 
latitude  stations.  The  theoretical  magnetograms  generally  agree  with  observations 
with  regard  to  the  largest  features,  although  the  model  is  far  from  being  sophisti¬ 
cated  enough  to  make  quantitatively  reliable  predictions.  Each  component  (,"3X, 

!\By,  ABZ)  of  each  theoretical  magnetogram  is  broken  down  into  individual  contribu¬ 
tions  from  ten  different  kinds  of  currents  (e.g.,  region-1  Birkeland  currents 
connecting  to  northern  hemisphere,  ring  currents).  All  major  magnetosphere- 
ionosphere  currents  are  now  included,  except  for  ground  currents  and  neutral-wind- 
priven  currents.  For  the  modeled  substorm  event,  the  change  in  vertical  component  i 
3Z  at  each  of  the  AFGL  stations  is  due  mainly  to  the  auroral  electrojet  and  is 
relatively  easy  to  predict  correctly.  The  change  in  the  eastward  component  By 
represents  the  small  net  effect  of  the  region-1  and  region-2  Birkeland  currents, 
which  nearly  cancel  each  other;  ABy  is  difficult  to  predict  with  quantitative 
accuracy.  The  change  in  northward  component  Bx  at  an  AFGL  station  during  a  sub¬ 
storm  represents  the  small  resultant  effect  of  many  different  kinds  of  currents 
(particularly  region-1  and  region-2  Birkeland  currents  from  both  hemispheres  and 
ring  current);  ABX  is  difficult  to  predict  with  quantitative  accuracy.  Our  general 
inclusions  concerning  interpretation  of  AFGL  magnetograms  in  substorms  are  the 
Following: 

(1)  The  standard  "substorm  current  loop"  (westward  electrojet,  reduced  tail  current!, 
:onnecting  Birkeland  current)  is  a  very  useful  mnemonic  for  remembering  the 
:haracteristic  ground-magnetic  disturbance,  but  it  does  not  really  represent  very 
well  the  physical  distribution  of  currents  flowing  during  a  substorm;  our  model 
(and  TRIAD  data)  suggests  that  the  actual  change  in  current  distribution  in  a  sub¬ 
storm  is  intensification  of  the  electrojet,  the  region-1  and  region-2  Birkeland 
currents,  and  the  ring  current, 

(2)  The  standard  picture  of  the  asymmetric  low-latitude  ground  disturbance  early 
in  a  magnetic  storm,  a  westward  duskside  partial  ring  current  completed  by  Birke¬ 
land  currents  and  an  eastward  electrojet,  also  is  a  useful  mnemonic  that  does  not 
represent  well  the  physical  currents  that  flow  during  ring-current  injection;  our 
nodel  suggests  that  the  actual  low-latitude  asymmetry  during  ring-current  injection 
is  due  primarily  to  strengthening  of  the  electrojets  and  of  the  region-1  and 
region-2  Birkeland  currents  that  feed  them,  and  not  to  a  drastic  dawn-dusk  asym¬ 
metry  in  the  ring  current; 

(3)  Our  model  results  suggest  that,  for  substorms  that  occur  at  times  such  that 
at  least  part  of  the  continental  U.S.  is  in  the  postmidnight  local-time  sector, 
the  ABz-components  of  the  magnetograms  from  the  northern-US  chain  of  AFGL  magneto¬ 
meters  might  provide  a  valuable  index  that  would  indicate  the  total  strength  of 
the  westward  auroral  electrojet,  particularly  if  combined  with  a  model  calculation 
that  corrects  for  ground  currents. 
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I.  INTRODUCTION 

In  the  early  1970' s,  it  became  clear  that,  in  standard  convection 
theory,  the  flow  of  plasma  and  electrical  currents  in  the  inner  magneto¬ 
sphere  and  corresponding  ionosphere  constituted  a  well-posed  theoretical 
problem;  this  problem  could  in  principle  be  solved  numerically,  given 
certain  nontrivial  boundary  conditions  and  input  parameters  (see,  e.g., 

Va syl iunas ,  1970,  1972).  There  has  been  a  longstanding  effort  at  Rice 
aimed  at  numerical  solution  of  this  inner-magnetosphere/ionosphere  con¬ 
vection  problem.  Over  the  years,  our  computer  modeling  activity  has 
gradually  improved  by  (i)  increasing  the  number  of  physical  processes 
included,  (ii)  imposing  more  realistic  inputs  and  boundary  conditions  and 
(iii)  performing  direct  comparisons  with  observations  (Wolf,  1970;  Jaggi 
and  Wolf,  1973;  Wol f ,  1974;  Hare!  and  Wolf,  1976;  Hare!  et  al ■ ,  1979; 

Hare!  et  al . ,  1980a,  b;  Spi ro  et  al . ,  1980).  By  the  mid-1970's,  it  was 
pretty  well  established  that  standard  convection  theory  was  at  least 
semi-quantitatively  consistent  with  the  basic  features  of  inner-magnetospheric 
currents,  electric  fields  and  plasma  flows.  For  example,  the  tendency  for 
low  L-values  to  be  shielded  from  the  convection  electric  field  and  the 
general  global  pattern  of  Birkeland  currents  were  correctly  predicted 
(Schield,  et  al ■ ,  1969;  Vasyliunas,  1972;  Jaggi  and  Wolf,  1973;  Wolf.  1974) 
before  they  were  observed  (e.g.,  Gurnett,  1972;  Hepjaner,  1972;  Zmuda  and 
Armstrong ,  1974). 

After  our  computer  model,  which  is  a  quantitative  embodiment  of  the 
standard  convection  picture,  proved  successful  in  predicting  these  major 
features  qualitatively,  it  seemed  ready  for  precise,  quantitative  confronta¬ 
tion  with  observations.  Since  magnetospheric  phenomena  are  highly  variable, 
precise  quantitative  comparison  with  observations  requires  theoretical 
computations  for  individual  events,  and  this  has  been  the  basic  thrust  of 
our  work  for  the  last  three  years. 

Our  work  under  contract  F19628-77-C-0005  has  had  two  major  parts: 

(1)  to  adopt  our  basic  computer  model  for  simulation  of  individual  events, 
using  real  event  data  as  input,  and  checking  model  results  against  other 
data  from  the  event(s);  (2)  to  do  the  extensive  programming  involved  in 
using  model  results  to  predict  individual  ground  magnetograms,  and,  by 
comparing  these  theoretical  magnetograms  for  the  chain  of  AFGL  magnetometer 
stations  with  the  actual  observations,  to  provide  a  unique  analysis  and 
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interpretation  of  those  observations. 

In  Section  II,  we  give  a  brief,  general  summary  of  the  work  we  did 
in  computer  simulating  events  in  the  inner  magnetosphere,  and  of  the 
results  and  physical  conclusions.  A  much  more  detailed  summary  is 
contained  in  three  long,  detailed  papers  describing  this  work,  submitted  as 
Scientific  Reports  under  contract  F19628-77-C-0005. 

In  Section  III,  we  summarize  the  work  we  did  in  computing  theoretical 
magnetograms,  give  a  complete  display  of  theoretical  magnetograms  for  the 
AFGL  chain  and  a  few  other  stations  (including  computed  contributions  to 
the  AB(t)  at  each  station  by  ten  different  kinds  of  magnetospheric  and 
ionospheric  currents),  and  our  detailed  analysis  and  interpretation  of 
the  comparison  with  observations. 

Section  IV  gives  a  chronological  list  of  the  various  programming  and 
data  analysis  tasks  that  were  carried  out  under  the  contract.  Section  V 
lists  contract-supported  reports  and  papers,  and  gives  infonnation  on 
fiscal  and  administrative  aspects  of  this  contracted  research  program. 
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II.  BRIEF  SUMMARY  OF  COMPUTER-SIMULATION  EFFORT 
A.  Model  Formulation  and  Logic 

Our  model  concentrates  on  the  inner  magnetosphere,  specifically  the 
region  where  magnetic  field  lines  are  certainly  closed  and  magnetic-f  old 
models  reasonably  reliable.  We  apply  our  boundary  conditions  at  L  '  10, 
which  maps  to  a  curve  slightly  equatorward  of  the  polar-cap-boundary,  and 
we  attempt  a  self-consistent  solution  equatorward  of  that  curve  (Fig.  1). 

Table  1  summarizes  the  assumptions  of  our  model  (specifically  the 
version  used  for  the  19  September  1976  event).  The  logic  of  the  same 
version  of  the  model  is  illustrated  in  Figure  2.  The  basic  logic  loop 
(the  central  pentagon  of  the  figure)  is  a  modification  of  a  diagram  given 
by  Vasyl iunas  (1970b).  Dashed  lines  mark  future  additions  to  our  program 
which  are  not  included  yet.  Starting  at  a  given  time  T  with  a  given  hot- 
particle  distribution  in  the  magnetosphere,  usually  estimated  from  average 
plasma-sheet  particle  data  (upper  box),  we  proceed  counter-clockwi se.  We 
first  compute  the  divergence  of  the  gradient-curvature  drift  current  in 
the  magnetosphere.  An  important  input  for  this  calculation  is  the  magnetic 
field  model.  For  the  19  September  1976  simulation  we  used  a  superposition 
of  the  01 son-Pfi tzer  (1974)  analytic  model  and  a  time-dependent  "substorm 
current  loop"  that  simulates  the  effects  of  an  induction  electric  field. 

From  the  divergence  of  gradient/curvature  drift  current  in  the 
equatorial  plane,  we  compute  the  divergence  of  ionospheric  current, 
assuming  that  Birkeland  currents  connect  the  ionosphere  and  equatorial  plane 
to  maintain  V*j  =  0.  Then  using  a  semi -empirical  model  of  height- integrated 
ionospheric  conductivity,  we  solve  the  current-conservation  equation  in  the 
ionosphere  for  the  ionospheric  potential  distribution  V.  Our  boundary 
conditions  are  the  following: 

(i)  Zero  electric  current  across  the  equator  (this  condition  follows  fr  hi 
the  assumed  symmetry  between  the  hemispheres,  a  reasonable  assumption 
for  19  September,  which  is  near  equinox;  actually  the  condition  we  apply 
is  that  of  zero  current  across  latitude  21’). 

(ii)  Specified  potential  V  on  the  "polar  boundary,"  the  equatorward  boundary 
of  the  rep i on- 1  Birkeland  current,  which  should  be  distinguished  from 
the  "polar-cap  boundary"  commonly  defined  by  electric  field  reversals. 
Because  of  the  i rregulati ties  of  the  polar  boundary  we  actually  speci- 
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fied  the  potential  on  a  circle  that  encompasses  the  polar  boundary. 
The  distribution  of  the  potential  around  this  curve  has  the  general 
form  suggested  by  Figure  1  of  Heel  is  et  al.  (1976);  the  magnitude 
of  the  potential  drop  was  estimated  from  real-time  S3-?  observations 
for  the  19  September  1976  event. 

Given  the  ionospheric  potential  distribution  we  use  the  magnetic- 
field  model  to  map  V  to  the  equatorial  plane.  In  the  present  simula¬ 
tion  we  neglect  the  component  E  that  is  parallel  to  B,  although  we  believe 
that  a  few  kV  potential  drop  over  limited  regions  would  not  affect  our 
results  substantially.  (One  of  our  priorities  for  the  near  future  is  to 
include  provision  for  parallel  electric  fields  in  our  model).  We  proceed 
with  the  logic  loop  to  calculate  magnetospheric  electric  fields.  The 
total  electric  field  is  the  sum  of  the  potential  electric  field  and  the 
induction  field  calculated  by  means  of  letting  the  equatorial  crossing 
point  of  the  field  line  vary  in  time. 

Given  the  potential  electric  field,  the  motion  of  equatorial  cross¬ 
ing  points  due  to  induction,  and  the  magnetic-field  model,  the  program 
calculates  total  drift  velocities  (E  x  B  +  gradient  +  curvature)  for 
plasma-sheet  particles.  Specifically,  it  computes  the  motion  of  the  inner 
edge  of  each  species  of  the  plasma-sheet  electron-ion  distribution. 

Loss  of  electrons  by  precipitation  was  included  in  the  19  September 
1976  model  by  making  a  conventional  assumption,  namely  that  the  electrons 
suffer  strong  pitch-angle  scattering.  Under  these  conditions,  the  inner 
edge  of  the  electron  plasma  sheet  is  often  essentially  a  precipitation 
boundary  jVasyliunas,  1968,  KenneU  1969). 

Given  the  veloci ties  of  different  components  of  the  inner  edge,  the 
program  advances  the  position  of  the  inner  edge  fir  each  component  by 
thr  amount  appropriate  to  one  time  step  t  (30  seconds  in  these  runs). 

This  closes  the  logical  loop  for  another  time  step,  and  so  on. 

In  the  actual  numerical  procedure,  the  program,  in  every  time  step, 
reinterpolates  the  magnetic  field  model,  recalculates  Birkeland  currents 
for  the  new  particle  and  magnetic-field  configuration,  reads  the  observed 
electron  fluxes,  readjusts  Pedersen  and  Hall  conductivities,  solves  the 
equation  for  ionospheric  potentials  using  a  21  x  28  grid,  reinterpolates 
the  mapping  to  the  equatorial  plane,  calculates  corotation,  curvature 


and  gradient  drifts,  recomputes  boundary  velocities  and  moves  the  inner 
edge  of  various  components  of  the  plasma-sheet  (which  we  represent  by 
400  boundary  points). 

The  problem  is  further  complicated  by  the  fact  that  the  inner  edge 
of  the  plasma  sheet  is  often  rather  thin  (of  the  order  of  one  grid 
spacing).  Electric  fields  can  vary  by  large  factors  through  this  edc^ 
region  and  often  change  sign.  In  other  words,  electric  fields  generated 
by  one  part  of  the  inner  edge  strongly  affect  particle  motions  in  other 
parts.  To  accurately  model  this  sub-grid-scal e  phenomenon,  we  have  had 
to  include  a  rather  intricate  sel f-correction  scheme,  which  substantially 
complicates  the  program. 

B.  Data  Comparisons 

Figure  3  shows  some  basic  observational  characteristics  of  the  sub¬ 
storm-type  event  of  19  September  1976,  and  Table  2  summarizes  results  of 
comparing  theoretical  predictions  for  the  event  with  observations.  The 
best  relevant  data  sets  for  this  event  were  those  from  the  Air  Force 
satellite  S3-2.  These  data  provided  the  primary  inputs  to  the  model 
and  by  far  the  most  illuminating  checks  of  theoretical  predictions. 

Ground  magnetograms,  particularly  those  from  the  AF6L  chain,  provided 
very  sensitive  and  detailed  checks  of  model  predictions,  which  are  dis¬ 
cussed  in  detail  in  Section  III.  Various  types  of  data  were  unavailable 
for  19  September  1976,  and,  in  several  such  cases,  we  compared  our 
theoretical  predictions  with  published  "substorm-average"  or  "typical- 
substorm"  data. 

In  addition  to  the  19  September  1976  substorm- type  event,  we  have 
developed  a  new  version  of  our  program  for  use  in  simulating  the  main 
phase  of  a  major  magnetic  storm,  such  as  the  ones  that  occurred  29-30 
July  1977  or  10-11  April  1978.  We  have  collected  available  data  for 
both  of  these  events.  We  did  do  a  simulation  run  through  the  first  few 
hours  of  the  29  July  1977  event.  We  still  have  some  numerica 1 -ana lysi s 
problems  with  this  new  magne t ic-s torm  version  of  the  program.  However, 
despite  these  numerical  problems  (which  affect  only  a  few  aspects  of  the 
model),  and  the  fact  that  we  only  ran  through  the  first  few  hours  of  the 
main  phase,  and  the  fact  that  we  have  been  able  to  make  only  a  few  compari¬ 
sons  with  data  so  far,  we  have  derived  some  significant  physical  results  al¬ 
ready.  We  used  our  theoretical -magnetogram  program  to  derive  a  theoretical  D 


index  as  a  function  of  time  through  the  event,  and  it  agrees  well  with 
observations;  this  means  that  our  first-try  model  of  the  event,  based 
on  auroral  conductivities  estimated  from  AE ,  and  on  a  magnetic-field 
model  compressed  as  indicated  by  solar-wind  data,  automatical ly  injected 
a  ring  current  with  a  strength  consistent  with  the  observed  Dst.  The 
computed  inner  edges  of  the  electron  and  ion  plasma  sheet  are  consistent 
with  data  from  one  pass  of  S3-3  (J .  F.  Fennel  1 ,  private  communication, 

1979),  but  computed  inner  edges  do  not  seem  to  agree  well  with  ATS-  6 
data  (T.  Moore,  private  communication,  1979). 

C.  Physical  Conclusions 

The  most  important  conclusion  of  our  computer-simulation  effort  to 
date  is  that  it  is  feasible  to  do  self-consistent,  time-dependent  modeling 
of  the  inner-magnetosphere/ionosphere.  It  has  also  proven  to  be  both  use¬ 
ful  and  practical  to  model  specific  events  and  to  compare  model  predictions 
with  observations.  Our  comprehensive  global  model  gives  substantial 
physical  insights  into  the  behavior  o'  the  magnetosphere  even  when  there 
is  disagreement  between  our  computed  results  and  specific  observations. 

A  few  of  the  physical  conclusions  we  have  reached  to  date  are  listed  below: 
(1)  There  is  a  strong  tendency  <-  •  the  region  equa forward  of  the  auroral 
zone  to  be  shielded  from  the  Mg  -  ,nt. itudr-  electric  field.  The  shielding 
was  partially  effective  durincj  most  of  our  modeled  substorm,  and  was  quite 
effective  in  the  recovery  phase.  In  our  simulation  of  the  29  July  1977 
storm,  however,  the  compression  of  the  magnetosphere  in  the  storm  sudden 
commencement  thoroughly  disrupted  the  shielding,  which  was  not  re-established 
for  more  than  an  hour.  This  disruption  of  the  shielding  allowed  injection 
of  the  storm-time  ring  current. 

(?)  The  substorm  simulation  predicted  the  one  observed  occurrence  of  rapid 
subauroral  plasma  flow,  supporting  an  analytic  theory  of  enhanced  subauroral 
flows  [Southwood  and  Wolf,  1978], 

(3)  Reasonable  ring  currents  were  predicted  by  our  model  for  both  the 
substorm  (19  September  1976)  and  magnetic  storm  (?9  July  1977)  cases. 
Conclusions  are  preliminary  because  only  short  time  periods  were  modeled 
(3-4  hours  after  onset).  However,  the  maximum  calculated  ring-current 
strengths  were  reasonably  consistent  with  observed  depressions  of  Pst  in 
each  case,  and  for  the  substorm,  the  normal  energy-dependent  ion-arrival 


times  £McIlwajn,  1974]  were  predicted  for  geosynchronous  orbit  near  dusk. 
Thus,  although  the  shielding  process  tends  to  inhibit  ring-current  in¬ 
jection,  model  electric  fields,  which  included  carefully  computed  shield¬ 
ing,  allowed  injection  of  apparently  realistic  ring  currents. 

(4)  Region-2  Birkeland  currents  are  generally  observed  to  extend  through¬ 
out  the  diffuse-auroral  region;  this  observational  fact  is  in  essential 
disagreement  with  our  model  region-2  Birkeland  currents. 

The  assumption  that  the  particles  gradient,  curvature  and 
E  x  B  drift  earthward  from  the  tail  implies  that  the  region-2  Birkeland 
currents  occur  at  the  inner  edge  of  the  plasma  sheet,  independent  of  most 
details  of  the  model.  This  discrepancy  may  well  be  pointing  out  a  signi¬ 
ficant  flaw  in  classical  convection  theory,  possibly  the  nonadiabatic  loss  of 
plasma  from  plasma-sheet  flux  tubes  [Erickson_  and  _Wol f ,  1 980]* 

(5)  The  total  Joule  heating  during  an  injection  event  is  comparable  in 
magnitude  to  the  change  in  ring-current  energy. 

(6)  The  simulation  results  have  proved  to  be  consistent  with  a  number  of 
heretofore  puzzling  established  observat ional  facts,  and  suggest  explana¬ 
tions  for  them.  These  include:  (i)  the  Birkeland-current  overlap  region 
near  midnight,  where  a  band  of  upward  currents  is  bounded  by  bands  of 
downward  currer ts  (observations:  I  i  jinia  and  Potemra ,  1978;  theory:  Harel 
et  al .  ,  1980b);  (ii)  the  shift  from  eastward  to  westward-directed  E  near 
midnight  at  L  =  4  in  substorms  (observations:  Carpenter  et  al . ,  1979; 
interpretation:  Spiro  et  aJL  ,  1980). 

Conclusions  based  on  theoretical  magnetograms  are  listed  separately 
in  Section  III. 


*  A  copy  of  the  paper  by  Erickson  and  Wo 1 f  (1980),  which  describes  work  that  was 
supported  in  part  by  contract  FI 9623-77-C-000S  but  was  not  submitted 
as  a  Scientific  Report,  is  included  as  an  Appendix. 
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III.  THEORETICAL  ANALYSIS  OF  GROUND  MAGNETOGRAMS 

Scientists  have  been  working  for  most  of  the  twentieth  century 
on  usinq  unique!  ic-f  ielci  measurements  t  rom  collections  of  qrourui 
stations  to  deduce  the  constantly  shifting  patterns  of  magnetospheric 
and  ionospheric  currents.  Rather  elaborate  mathematical  models  of 
"equivalent  ionospheric  currents"  were  constructed.  This  type  of  re¬ 
search  progressed  much  more  rapidly  in  the  sixties  and  seventies,  as 
more  and  more  rocket  and  spacecraft  data  could  be  compared  with  the 
g round  data,  eliminating  some  of  the  ambiguity  resulting  from  the  fact 
that  the  field-aligned  Birkeland  currents,  which  drive  most  horizontal 
ionospheric  currents,  are  hard  to  measure  by  means  of  ground  magneto¬ 
meters.  Our  approach  to  the  analysis  of  ground  niagnetograms ,  under 
contract  F19628-77-C-0005,  has  been  to  use  our  computer  simulation  of  a 
substorm  event,  based  on  some  satellite-gathered  input  data  and  a  great 
deal  of  theory  and  computation.  We  compute  a  reasonably  realistic, 
complicated,  3-dimensional  current  system,  and  then  compute  ground- 
magnetic  variations  by  straightforwardly  integrating  the  Biot-Savart  law 
over  this  current  system.  The  results  are  then  compared  with  the  actual 
data,  and  the  results  analyzed. 

A.  Theoretical  Computation  of  Ground  Magnetograms 
Running  our  main  simulation  program  through  an  event  furnishes  a  time- 
history,  for  the  event,  of  the  magnetospheric  plasma  distribution  within  our 
modeling  region  and  of  ionospheric  and  maqnetospheric  electric  fields. 
Time-dependent  magnetic-field  and  ionospheric-conductivity  models  are 
used  as  input.  Using  all  of  these  model  outputs,  permanently  stored 
after  the  main  simulation  run,  we  compute  the  following  currents: 

(i)  north-south  ionospheric  currents  between  grid  points,  assumed  equally 
divided  between  the  two  hemispheres; 

(ii)  east-west  ionospheric  currents  between  grid  points,  also  assumed  equally 
divided  between  the  two  hemispheres; 

(iii)  region-2  Birkeland  currents  connecting  ionospheric  grid  points  to  grid 
points  in  the  equatorial  plane;  these  currents,  divided  equally 
between  the  two  hemispheres,  are  computed  by  satisfying  Kircbhoff's 
law  at  ionospheric  grid  points; 
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(iv)  ring  currents  out  in  the  magnetosphere,  computed  from  the  plasma 
pressure  gradients  indicated  by  the  simulation;  for  simplicity,  the  total 
ring  current,  integrated  along  magnetospheric  flux  tubes,  from  the 
southern  ionosphere  to  the  northern  ionosphere,  is  represented  as  flowing 
along  wires  in  the  equatorial  plane. 

To  obtain  a  reasonable  representation  of  the  global  current  system,  wn 

include  the  currents  in  the  region  just  poleward  of  our  high-latitude 

boundary,  which  is  defined  to  be  the  eguatorward  boundary  of  the  region-1 

currents.  Using  measured  S3 -2  electron  fluxes,  we  construct  a  band,  poleward 

of  our  main  simulation  boundary,  representing  the  region  of  aurorally  enhanced 

conductivities,  where  large  horizontal  currents  flow.  The  band  is  taken  to 

have  latitude-independent,  local-time  dependent  Hall  and  Pedersen 

conductivities  [£..  and  l  ].  The  values  of  L'  and  band  width  on  the 
v  H  pJ  H  p 

subsatellite  track  are  chosen  as  a  best  fit  to  measured  S3-2  electron 
fluxes.  The  current  into  the  equatorward  boundary  of  this  high-latitude 
auroral  band  at  any  local  time,  is  taken  to  be  equal  to  the  value  computed  by 
the  main  simulation  program.  The  current  into  the  poleward  boundary  of  the 
high-latitude  auroral  band  is  taken  to  be  zero.  That  is,  the  polar  cap  is 
taken  to  be  an  insulator.  We  assume  that  region-1  Rirkeland  currents  are 
independent  of  latitude  within  the  band. 

Figure  4  shows  qualitative  sketches  of  our  total  three-dimensional  model 
current  system,  which  involves  about  2700  wires  or  bands.  We  integrate  the 
Biot-Savart  law  numerically  over  this  maze  of  wires,  at  20-minute  intervals 
through  the  modeled  event.  Actually,  for  wires  near  the  ground-observing 
point,  we  replace  the  wires  by  appropriate  bands,  and  assume  that  the  current 
is  uniformly  distributed  over  the  band,  which  is  centered  on  the  original 
wire.  North-south  horizontal  bands  have  widths  equal  to  the  spacing  of  grid 
points  in  local  time,  fast-west  horizontal  bands  have  widths  equal  to  local 
latitudinal  grid  spacing.  Birkeland-current  bands  are  oriented  roughly 
perpendicular  to  the  local  meridional  plane  and  have  widths  equal  to  the 
spacing  of  grid  points  in  local  time. 
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R.  Numerical  Results  and  Discussion 

Wo  present  in  Figures  5-25,  detailed  numerical  results  from  our 
thooretical-magnetograni  computations,  there  are  usually  four  figures  for  each 
station:  the  first  three  give  the  contributions,  to  each  component 
of  All,  from  each  of  ten  types  of  current  ;  the  fourth  shows  the  comparison 
between  theory  and  observations. 

Por  the  AFGL  stations,  three  curves  are  shown  for  each  comparison  between 
theory  and  observat ions :  (1)  a  solid  curve,  representing  the  theoretical 

prediction  for  AB  due  to  magnetosphere-dri ven  currents;  (2)  a  dashed  curve, 
representing  data  for  19  September  1970;  (3)  a  dotted  curve  representing  our 
estimate  of  what  AR  would  have  heen  on  19  September  1976,  if  there  had  been  no 
magnetospheric  activity.  The  dotted  curve  is  based  on  several  magnetograms 
from  the  same  station,  from  very  quiet  days  in  late  summer  and  early  fall. 

The  theoretical  curve,  representing  magnetospheric  effects  should  in  principle 
agree  with  the  difference  between  the  dashed  and  dotted  curves. 

Newport ,  Washington  Station,  x-component  (northward)  (figure  6;  hot t on 
panel  of  figure  8).  There  is  poor  correlation  between  theory  and  data  here. 
The  main  reason  for  the  difficulty  of  making  accurate  computations  is  easily 
seen  in  figure  6.  Many  kinds  of  currents  contribute  appreciably  to  AB^  at 
Newport.  For  example,  the  types  of  current  that  typically  make  contributions 
exceeding  one-third  of  the  total  resultant  are  the  following:  region-1 
Birkeland  curronts-northern  and  southern  hemispheres  (Rl),  region-2  Rirkeland 
currents-northern  and  southern  hemispheres  (R ?),  and  magnetospheric  ring 
current  (RC).  The  sun  of  the  absolute  values  of  the  various  contributions  is 
typically  5-10  times  the  actual  total.  Since,  in  our  type  of  model,  factor- 
of-two  errors  in  individual  currents  could  easily  occur  due  lo  inadequate 
input  data,  there  is  little  hope  of  our  making  reliable  predictions  of  .Ml  . 
I'ven  from  a  non-modeling  point,  of  view,  the  AB  measured  at  Newport  is  hard 
to  interpret.,  since  it  is  such  a  complicated  combination  of  different  kinds  of 
currents. 

Newport  y-componont.  (eastward)  (T 1  gore  6j_  midd  I e  panel  of  F  igure  B) .  I n 
this  case,  there  is  a  resemblance  between  the  theoretical  curve  and  the  real 
one.  However,  the  model  failed  to  reproduce  the  two  observed  large  dips 
centered  at  about  1055  UT  and  1220  '.IT;  t  he  model  oroaMy  underest  i  mated  t  ho 
amplitude  of  the  first,  dip,  and  cornel  et  >ly  missed  t  he  second  one.  The 
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computed  AB  component  is  essentially  the  result  of  competition  between  the 
oppositely  directed  region-1  and  region-2  currents.  The  AB^  tends  to  be 
positive  in  the  event,  indicating  that  the  region-2  currents  generally  fend  to 
have  larger  effect,  because  they  are  nearer  to  the  station.  There  are  two 
problems  with  making  accurate  model  predictions  of  AB  .  First,  there  is 
relatively  close  competition  between  opositely  directed  region-1  and  region-2 
currents,  which  makes  accurate  prediction  of  the  sum  difficult.  The  maximum 
discrepancy  between  computed  and  observed  total  AB  's  is  about  30y,  which  is 

y 

not  terribly  large,  considering  that  the  AB  's  from  each  type  of  Birkeland 

current  average  about  80y  absolute  magnitude.  Second,  the  nature  of  our 

basic  input  data,  which  come  from  polar  passes  of  a  satellite  and  therefore 

have  an  effective  time  resolution  of  about  50  minutes,  prevents  accurate 

simulation  to  the  time  development  of  a  substurm.  Judging  from  this  Newport 

magnetogram  and  several  magnetograms  from  stations  further  west,  there  was  a 

brief,  localized  substorm  with  peak  about  12  UT,  centered  at  Paci f ic-Ocean 

longitudes.  Our  low-time-resolution  input  data  caused  us  to  combine  that 

substorin  with  the  recovery  phase  of  the  first  substorm.  We  also  note  that 

positive  AB  developed  in  the  hour  before  onset  of  the  10  UT  substorm--to  a 
y 

much  greater  extent  in  the  data  than  in  the  nodel.  This  may  indicate  an 
under-estimate  of  the  strength  of  the  growth  ohase.  It  is  also  interesting  to 
compare  the  signature  that  would  be  predicted  for  AB  using  the  standard 
substorm  current  loop  (e.g.,  McPherron,1973) ,  which  involves  downward 
Birkeland  on  the  dawnside  during  the  substorms.  That  picture  would  predict 
the  negative  excursion  in  By  near  the  peaks  of  the  two  substorms,  but  would 
not  predict  the  tendency  for  By  with  19  September  event  to  be  greater,  as  the 
average,  than  on  a  quiet  day.  In  our  model  the  positive  By  was  due  to 
strengthening  of  region-2  currents  during  the  event,  in  any  case,  the  simple 
substorm  current  loop  gives  a  misleadingly  simple  picture  of  Birkeland 
currents  flowing  in  a  substorm.  Actually,  there  are  two  oppositely  directed 
sets  of  Birkeland  currents,  each  of  which  varies  substantially  during  the 
substorm. 

Newport  /-component  (downward)  (Figure  7 ;  top  panel  of  Figure  B) . 

The  AB  vector  at  Newport  was  basically  upward,  both  in  the  data  and  the  model, 
and  the  model  predicted  this  dominant  part  of  AB  with  remarkable  accuracy. 

The  difference  between  theory  and  observation  in  this  case  is  less  than  the 
effect  of  ground  currents,  which  aren't  included  in  the  model  and  must  be 
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appreciable.  It  is  clear  from  Figure  7  that  many  kinds  of  currents  contribute 

significantly  to  ABz  at  Newport,  particularly  region-1  and  region-2  Birkeland 

currents  (R1  and  R2),  north-south  ionospheric  current  (N-S),  ring  current 

(RC),  east-west  ionospheric  current  inside  our  modeling  region  (F-W)  and  east- 

west  electrojet  current  poleward  of  our  modeling  boundary  (EJ).  Southern- 

hemisphere  currents  fortunately  do  not  contribute  much.  Fukushima's  theorem 

( Fukushiina ,  1969;  Vasyl i unas ,  1970a]  states  essentially  that  ground  magnetic 

effects  of  Birkeland  currents  and  Pedersen  currents  cancel  each  other  out, 

under  various  simplifying  assumptions,  particularly  uniform  Pedersen  currents 

and  straight,  vertical  Birkeland  currents.  If  one  ignores  inaccuracies 

involved  in  these  assumptions,  which  are  not  made  in  our  computer  simulations, 

one  would  expect  approximate  cancellation  of  the  magnetic  effects  of  the 

region-1  and  region-2  Birkeland  currents,  and  the  Pedersen  parts  of  the 

horizontal  ionosphere  currents.  In  our  model,  most  of  the  north-south 

ionospheric  current,  part  of  the  F-W  current  in  our  modeling  region,  and  a 

significant  part  of  the  electrojet  poleward  of  our  modeling  boundary,  are 

Pedersen  currents.  Naive  application  of  Fukushima's  theorem  to  our  case  would 

suggest  that  the  ground  magnetic  variation  is  due  to  the  majority  Hall-current 

part  of  our  east-west  ionospheric  current  (F.W  +  FJ),  plus  the  ring  current 

(RC).  The  actual  computed  total  AB  in  Figure  7  is  approximately  equal  to 

RC  +  0./5  (FJ  4-  FW).  The  reason  that  AB?  is  easier  to  predict 

than  AB  or  AB  is  that  the  Hall  current,  which  does  not  tend  to  cancel  out 
x  y 

according  to  Fukushima's  theorem,  gives  a  large  AB  at  65°  geomagnetic 
latitude  at  night. 

Our  model  results  suggest  that  AB  measurements  from  a  chain  of  stations 
near  bb°  geomagnetic  latitude  might  furnish  a  useful  estimate  of  total 
electrojet.  strength  --  more  global  than  the  AF  index,  which  is  sensitive  to 
conditions  right  over  the  stations  used.  Stations  near  bb’,  being  further 
from  the  electro  jet,  night  furnish  a  better  estimate  of  total  strength.  Of 
course,  to  obtain  quantitatively  accurate  information  on  electrojet  strength, 
corrections  would  have  to  be  made  for  the  ring-current  effect  (which  can  be 
independently  es* mated  from  Dst)  and  for  ground  currents. 

Camp  Douglas  x-component .  (I i gure  9 ,  bottom  panel  of  Figure  12). 

Same  comments  as  for  Newport,  x-componer  t. . 

Cam;)  Douglas  y-conponent .  (F i gurt  10,  middle  panel  of  F i gure  12). 


Same  comments  apply  as  for  the  newport  y-magnetogram,  except  that  in  this  case 
the  substorm  at  **  12  UT  did  not  show  up  clearly  in  the  observations, 
presumably  because  the  station  was  too  far  east.  Also,  in  this  case  the 
positive  average  value  of  AB  is  probably  due  to  quiet-time  wind-driven 
currents,  not  to  magnetospheric  effects. 

Camp  Oougl as  z-coinponent .  (Figure  11 ;  top  panel  of  Figure  12) . 

Same  comments  apply  as  for  the  Newport  z-magnetogram.  The  discrepancy  in 
magnitudes  of  AR^  is  well  within  expected  modeling  errors.  The 
computed  A8  is  approximately  given  by  RC  +  CJ  ^  [W. 

Sudbury  x-component.  (Figure  13;  bottom  panel  of  Figure  16) . 

Agreement  is  adequate  in  this  case,  hut  of  marginal  significance,  since  the 
net  A8^  is  small  and  made  up  of  many  comparable  contributions. 

Sudbury  y-coinponents  (Figure  14 ;  middle  panel  of  Figure  16) .  Agreement 
between  theory  and  observations  is  poor,  and  the  observed  magnetogram  is 
difficult  to  interpret,  except  for  the  substorm-associated  minimum  just  before 
1100  UT  .  Quiet-day,  wind-driven  currents  are  relatively  large,  because  the 
station  is  in  sunlight  for  most  of  the  event. 

Sudbury  z-component  (Figure  15;  top  panel  of  Figure  16).  Agreement  is 
again  very  good.  The  relation  ABz=  RC  +  (CJ  +  FW)  works  well  early  in  the 
event.  Later  in  the  event,  when  the  station  is  well  past  dawn,  effective 
electrojet  strength  drops  and  deviations  from  Fukushima's  theorem  become 
relatively  more  important. 

Theoretical  Magnetograms  for  "Equatorial  Stations" 

Figures  17-20  show  what  would  be  seen  at  four  imaginary  "stations"  that 
sit  on  the  geomagnetic  equator  at  constant  magnetic  local  time.  We  plot  only 
x-components,  because  y-  and  z-components  vanish  by  symmetry  in  our  model. 
Equatori al -electrojet  currents  are  not  included  in  our  model,  so  these  results 
are  more  representati ve  of  what  would  be  observed  off  the  equator,  but  at 
la!  i tudes  <  3(1° . 

In  every  case,  the  ring  current  makes  the  largest  contribution 
to  |AB  |.  However,  on  the  dawn  side,  the  combination  of  east-west, 
ionospheric  currents  and  a  partial  non-cancel lation  of  Birkeland  currents  acts 
to  increase  B  .  Altogether,  the  depression  of  Bx  is  less  than  half  as  great 
on  the  dawn  side  as  on  the  dusk  side.  The  main  reason  for  the  dawn-dusk 
asymmetry  in  AR  is  that,  ’n  the  simulations,  there  is  a  net  downward 
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Hi rkel and  current  on  the  day  side  and  a  net  upward  current  on  the  night  side 
(Hare]  et  al  . ,  1980b).  Correspondingly,  there  is  a  net  current  flow  in  the 
electrojets  from  day  side  to  night  side.  Thus  there  is  an  effective  high- 
latitude  current  from  day  side  to  night  side  in  each  hemisphere,  causing  a 
depression  of  Bx  on  the  dusk  side,  enhancement  on  the  dawn  side.  This,  we 
propose,  is  the  main  cause  of  the  low-latitude  dawn-dusk  asymmetry  in  AB  near 
the  equator.  The  net  Bi rkel and  current  down  on  the  dayside,  up  on  the  night 
side,  should  not  be  visualized  as  being  completed  by  a  simple  westward  partial 
ring  current  on  the  dusk  side.  The  model  ring  current  is,  in  total  strength, 
nearly  symmetric  about  local  midnight,  although  some  components  (e.g.,  30  keV 
ions)  may  be  stronger  at  dusk  than  at  dawn  (or  vice  versa).  Region-1 
Birkeland  currents  are  completed  by  a  complicated  and  unknown  combination  of 
tail  current,  Chapman-Ferraro  currents  and  megnetosheath  currents,  while 
region-2  Birkeland  currents  are  completed  mostly  by  gradient-and  curvature- 
drift  currents  on  closed  field  lines.  It  is  misleading  to  think  of  the  net 
downward  Birkeland  current  on  the  day  side  (an  incomplete  cancellation  of 
region-1  and  region-2  currents)  and  corresponding  net  upward  current  on  the 
night  side,  as  being  linked  by  a  simple  partial  ring. 

Figure  21  compares  the  average  total  aB^  from  Figures  17-20  with  the 
observed  Dst.  In  this  local-time  averaging  process  for  equatorial  stations, 
the  contributions  from  Birkeland  currents  and  horizontal  ionospheric  currents 
cancel  out,  leaving  just  the  ring-current  effect.  For  the  first  hour  or  two 
after  substorm  onset,  the  ring  current  injected  in  our  simulation  depressed 
the  equatorial  field  by  about  the  amount  observed  in  Dst.  However,  our  ring 
current  energy  continued  to  increase  (Hare!  et  al.,  1980b)  and  the 
"theoretical  Dst"  shown  in  Figure  21  continued  to  decrease  until  1300.  This 
disagrees  with  the  observations,  which  showed  Dst  starting  to  recover  between 
1200  and  1300  UT.  The  observed  average  Dst  for  the  hour  1300-1400  UT  was  near 
its  initial  value,  indicating  a  rapid  ring-current  recovery  that  our  model 
will  probably  have  difficulty  reproducing  if  the  program  is  run  another  hour. 

We  present  one  example  of  an  auroral-zone  station,  namely  Fort  Churchill, 
Canada . 

tort  Churchi  1 1  x- component .  (Figure  22 ;  bottom  panal  of  Figure  2b) . 
Agreement  is  very  acceptable.  The  mod>‘l  overest.  imated  the  maximum  depression 
of  Bx,  but  by  a  factor  that  is  well  wi  hin  expected  modeling  errors. 
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Fort  Churchi 1 1  y-component.  (Figure  23;  middle  panel  of  F i gure  25) . 

There  was  little  meaningful  y-perturbation  in  either  the  theory  or  the  data. 

As  expected  on  the  basis  of  Fukushima's  theorem,  contributions  due  to  region 
1,  region  2,  and  north-south  ionospheric  currents  almost  cancel. 

Fort  Churchi 1 1  z-component.  ( F i gure  24;  top  panel  of  Figure  25) . 
Agreement  here  is  again  very  good.  The  downward  field  is  due  to  the  station 
being  north  of  most  of  the  westward  electrojet. 

Theoretical  magnetograms  have  been  computed  for  two  other  stations.  The 
theory  and  data  for  Boulder,  Colorado  are  too  similar  to  Newport  and  Camp 
Douglas  to  be  independently  interesting.  The  theory  and  data  disagree  for  the 
College,  Alaska  station,  due  to  a  known  problem  with  the  model:  during  the 
main  part  of  the  substorm.  College  lies  in  the  model's  Harang-di sconti nui ty 
region,  which  was  too  far  east  due  to  a  poor  choice  for  the  potential  on  our 
poleward  boundary.  The  problem  had  earlier  been  identified  in  comparisons 
with  TRIAD  data. 


C.  Conclusions  from  Theoretical  Magnetograins  and  from 
Comparisons  With  Data 

Our  major  conclusions  from  this  generation  of  theoretical  magnetograins 
and  comparison  of  them  with  data  are  the  following: 

(L)  The  agreement  between  initial  simulation  results  and  observations  have 
turned  out  to  be  good  enough  that  one  can  learn  a  lot  from  the  comparison. 

(2)  However,  to  attain  good  agreement  between  our  simulations  and  substorm¬ 
time  magnetograms  will  require  at  least  good  input  at  <  20-minute  intervals, 
which  cannot  be  achieved  using  a  single  polar-orbiting  satellite.  The 
difficulty  of  obtaining  appropriate  high-time-resolution  input  data  is  one  of 
our  fundamental  problems. 

(3)  Our  simulation-based  theoretical  magnetograms,  which  involve  fairly 
realistic  3-dimensional  current  systems  and  detail  the  contributions  of 
different  kinds  of  currents  to  the  AB's  at  various  ground  stations,  provide  a 
now  appreciation  for  the  complexity  of  the  source-current  distributions  for 
observed  ground-magnetic  variations. 

(4)  Although  the  standard  "substorm  current  loop"  (disruption  of  tail 
current,  westward  electrojets,  connecting  Rirkeland  current)  represents  part 
of  the  change  in  current  distribution  in  the  expansion  phase  of  a  substorm,  it 
is  more  accurate  to  describe  the  substorm  current  as  enhancement  in  the  entire 
region- 1  and  region-2  current  systems,  and  their  associated  Pedersen  and  Hall 
currents.  The  substorm  current  loop  does  not  necessarily  represent  the 
physically  dominant  current  changes  in  a  substorm,  only  the  ones  that  are  most 
obvious  from  the  ground. 

(5)  Our  computed  local -time-averaged  low-latitude  depression  of  B  agreed  well 
with  the  observed  Dst  for  the  first  hour  or  so  after  substorm  onset.  However, 
our  theoretical  Dst  continued  to  indicate  a  slow  increase  in  ring-current 
strength  after  that,  while  the  observations  indicated  decreasing  strength.  We 
do  not  have  an  explanation  for  this  discrepancy. 

(o)  The  observed  dawn-dusk  asymmetry  in  the  depression  of  Bx  at  low  latitudes 
during  ring  current  injection  is  not  due  primarily  to  a  physical  partial  ring 
current  injected  onto  the  dusk  side  but.  rather  to  a  net  downward  Birkeland 
current  on  the  dayside,  net  upward  current,  on  the  night  side,  and  generally 
ant i sunward  flowing  electrojets.  These  net  downward  dayside  and  upward 
nightside  Birkeland  currents  are  complicated  combinations  of  region-1  and 
region-2. 
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(7)  The  x-components  of  substorm  magnetograms  near  55°  latitude  are  due  to  a 
complicated  combination  of  currents,  and  cannot  be  accurately  modeled  yet. 

(8)  The  y-components  of  substorm  magnetograms  hear  55°  latitude  are  due  to 
imbalance  of  region-1  and  region-2  Birkeland  currents.  Since  AB^  is  due  to  a 
small  difference  between  two  large  quantities,  it  is  hard  to  predict  with 
quantitative  accuracy. 

(9)  On  the  night  side  the  2-components  of  55°  substorm  magnetograms  are 
basically  indicative  of  electrojet  strength.  The  model  was  successful  in 
predicting  these.  A  chain  of  stations  near  58"  geomagnetic  latitude  might 
provide  a  useful  index  of  electrojet  strength. 


IV.  Recapitulation  of  the  Work  Fffort 


This  section  recapi tulates ,  in  approximate  chronological  order,  the 
individual  steps  in  the  contracted  research,  as  mentioned  in  the  quarterly 
reports.  The  items  mentioned  below  mainly  concern  programming  and  data 
acquisition  and  display.  The  various  reports  on  the  work,  for  the  scienti¬ 
fic  community,  are  listed  separately  in  Section  V.  Scientific  conclusions  and 
results  are  summarized  in  sections  II  and  III. 

1.  A  Calcomp  plotting  routine  was  developed  for  automatic  plotting  of 
equipotential  patterns  and  particle  boundaries,  in  the  ionosphere  and 
equatorial  plane. 

2.  The  main  program  was  run  for  more  than  seven  hours  magnetosphere  time, 
mainly  as  part  of  the  procedure  of  testing  for  programming  errors,  for 
numerical  accuracy,  and  to  maximize  speed  and  minimize  core  size. 

3.  Data  were  acquired  from  the  Air  force  and  from  the  World  Data  Center 
for  use  in  choosing  an  event  to  simulate. 

4.  Programming  was  completed  for  inclusion  of  electron  precipitation  in  the 
program. 

5.  Programming  was  completed  for  computation  of  the  magnetic  field  in  a 
self-consistent  (through  approximate)  way  that  includes  pressure  balance 
at  the  inner  edge  of  the  plasma  sheet.  Inis  self-consistent  magnetic- 

field  routine  was  abandoned  eventually:  the  results  for  the  so  1 f-cons i stently 
computed  magnetic  field  were  not  greatly  di:  erent  from  results  with  the 
01 son-Pfi tzer  model  ,  but  the  routine  generated  bothersome  numerical  noise. 

6.  The  19  September  1976  event  was  selected  tm  simulation. 

7.  Data  from  the  S3-2  electric-field  detector  ae  1  electron  detector,  were 
acquired  tor  the  19  September  1976  event,  along  with  data  from  the  AFGL 
Magnetometer  chain  and  a  few  other  sources. 

b.  The  S3-2  electric -field  data  were  used  to  estimate  the  polar-cap  potential 
drop  as  a  function  of  time  through  the  event.  Measured  electron  fluxes  were 
used  to  compute  a  global,  time-dependent  conductivity  model. 

9.  The  substorm  current  loop  (westward  electro. let  ,  reduction  of  tail 
current,  connecting  current)  was  included  in  the  magnetic-field  model  used 
as  input  for  our  main  simulation  program. 

10.  The  first  computer-simulation  run  through  the1  event  of  19  September  10/6 
('vent,  was  completed.  Computed  electric  fields  were  compared  with  S3-!1  data. 

11.  Programming  was  completed  for  plotting  of  fine-structure  electric  fields. 


for  testing  of  numerical  accuracy  and  more  accurate  information  on  model 
electric  fields. 

12.  The  original  simulation  was  rerun  using  a  slightly  modified  numerical 
method  for  computing  the  fine-structure  electric  fields  and  a  smoother 
conductivity  profile. 

13.  The  19  September  1976  simulation  was  rerun  without  the  "substorm 
current  loop"  included  in  the  input  magnetic  field  model. 

14.  The  19  September  1976  simulation  was  rerun  with  a  higher  maximum  assumed 
potential  drop  (145  kV). 

15.  Programming  was  completed  for  following  the  motion  of  an  arbitrary 
drifting  particle,  using  the  main-program  results  for  electric  fields  as 
functions  of  time  through  the  event. 

17.  Programming  was  completed  and  tested  for  automatic  plotting  of  Birke- 
land  current  strengths,  as  functions  of  latitude  for  given  local  time. 

The  results  were  compared  with  average  TRIAD  data. 

18.  Programming  was  completed  for  making  arrow  plots  of  electric  fields 
and  flow  velocities  in  the  ionosphere  or  equatorial  plane. 

19.  Ring-current  energy  and  Joule  heating  were  computed  and  compared  for 
the  19  September  1976  event. 

20.  Programming  was  completed  for  computing  transverse  AB‘s  caused  by 
Birkeland  currents  along  the  S3-2  orbit.  Results  were  compared  with  actual 
S3-2  magnetometer  data. 

21.  A  large  program  was  written  and  tested  for  computing  ground-magnetic 
variations  by  integrating  the  Biot-Savart  law  over  a  maze  of  wires  represent¬ 
ing  the  currents  computed  in  our  main  simulation. 

22.  An  initial  set  of  theoretical  magnetograms  was  computed,  for  the  AFGL 
chain  and  several  other  magnetometer  stations.  Results  were  compared  with 
data. 

23.  Electric  fields  at  L  =  3-4  were  computed  and  compared  with  data  from 
whistlers  and  incoherent-backscatter  radars. 

24.  Analytic  models  were  constructed  for  the  high-latitude  portion  of  the 
auroral  zone,  poleward  of  the  boundary  of  our  main  simulation.  Results 
were  converted  into  a  form  suitable  for  inclusion  in  the  program  for 
computing  ground  magnetic  variations. 

25.  The  speed  and  accuracy  of  our  program  for  computing  ground  AB '  were 


increased  by  approximating  the  current  carriers  near  the  observation  point 
as  bands  rather  than  wires. 

26.  Analytic  work  was  carried  out  to  treat  the  pressure  variations  in 
converting  plasma -sheet  flux  tubes  out  in  the  tail.  (A  copy  of  the 
draft  paper  on  this  topic  is  reproduced  in  the  Appendix.) 

21 .  Data  were  acquired  for  the  29  July  1977  and  10  April  1978  magnetic 
storms . 

28.  Our  main  program  was  converted  into  a  form  suitable  for  a  magnetic 
storm  (including  a  magnetospheric  compression  and  formation  of  a  full 
ring  of  current) . 

29.  An  initial  computer-simulation  run  was  made  through  the  29  July  1977, 
and  a  few  dati  comparisons  were  made. 

30.  A  new  generation  of  theoretical  magnetograms  was  computed  for  the 
AFGL  chain  and  several  other  station,  ior  19  September  1976.  These  new 
magnetograms  integrate  over  current-carrying  bands,  rather  than  wires, 
near  the  observing  station,  and  include  ionospheric  and  Birkeland  currents 
poleward  of  the  boundary  of  all  main  simulation  program,  results  were 
compared  with  data  from  the  AFGL  chain  and  other  stations. 

Various  facets  of  this  computer-simulation  effort  have  been 
supported  by  several  grants  and  contracts  over  the  last  3]  years:  NSF 
grants  ATM74-21185  and  ATM79-2017,  which  have  supported  the  majority 
the  basic  program-development  effort  and  also  supported  the  simulation 
of  the  29  July  1977  event;  NASA  grant  NGR-44-006-137 ,  which  supported  a 
substantial  part  of  the  data  analysis;  Air  Force  contract  F19628-77-C-0012 , 
which  supported  the  early  conductivity  and  plasmasphere  modeling;  and 
Air  Force  contract  FI 9628-78-C-0078,  which  supported  analysis  of 
S3-2  data.  Contract  F19628-77-C-0005  supported  all  work  on  theoretical  and 
observed  magnetograms  done  from  1977  until  September  1979,  some  work  on 
S3-2  data  analysis,  and  part  of  the  basic  program-development  work  required 
to  carry  out  the  analysis  and  i nterpretation  of  Air  Force  data.  The  final 
part  of  the  analysis  of  AFGL  magnetograms  was  supported  by  Rice  University. 
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Abstract.  We  present  a  theoretical  argument  suggesting  that  steady, 
adiabatic  convection  probably  cannot  occur  throughout  a  closed-magnet ic-f iel.l- 
line  region  that  extends  into  a  long  magnetotail.  The  argument  is  applied 

quantitatively  to  the  earth's  magnetosphere,  using  several  magnetic-field 
models,  most  based  on  averaged  observations.  We  find  that,  if  there  were 
slow,  steady,  sunward,  adiabatic  convection  across  most  of  the  width  :  f  the 
plasma  sheet,  particle  pressure  would  increase  much  too  fast  with  decreasing 
geocentric  distance  to  be  consistent  with  tail -lobe  observations.  We 

hypothesize  that  sunward  convection  must  necessarily  be  time-dependent,  and 

that  the  magnetospheric  substorm  constitutes  the  essential  time-dependen* 
process  in  which  plasma-sheet  flux  tubes  lose  plasma  so  that  they  can  convoct 
sunward  to  the  i nner-pl asma-sheef  region  and  eventually  to  the  dayside 

magnetosphere. 


i 


Introdu'  tion 


The  solar  wind  and  the  earth's  magnetosphere  both  vary  greatly  with  time, 
and  it  is  clear  that  many  changes  in  the  solar  wind  trigger  changes  in  the 
magnetosphere.  On  the  other  hand,  it  is  certainly  not  clear  that  all 
substantial  magnetospheric  changes  are  caused  by  specific  solar-wind 
changes.  Workers  in  the  field  have  often  asked  themselves  (and  each  other! 
the  following  hypothetical  questions:  "it  the  solar  wind  incident,  on  the 

earth's  magnetosphere  were  held  completely  steady,  would  the  magnetosphere 
still  exhibit  gross  .  time  variations?  Specifically,  would  there  still  be 
magnetospheric  substorms?"  One  purpose  of  this  paper  is  to  argue 
theoretically  that  the  answer  to  both  of  these  questions  is  "yes." 

Cowley  (19/8)  and  Schindler  ( 1 97'?)  have  discussed  the  difficulty  of 
finding  se 1 1 -cons i stent  solutions  to  the  problem  of  slow,  t imo- i ndependont 
plasma  convection  throughout  a  closed-field-line  region  that  extends  well  out 
into  ,i  long  magnetotail.  We  argue  that  such  solutions  probably  do  not 
exist,  lor  the  case  of  the  earth  and  the  isotropic-pressure  uppruximat ion ,  we 
show  that  the  idea  of  slow,  steady,  uni  form,  sunward  convection  in  the  earth's 
plasma  sheet  is  clearly  inconsistent  wi'n  certain  magnet ic-1 ie I d  models,  most 
based  on  averaged  observations.  This  result,  suggests  a  mild  dilemma,  since 
there  is  overwhelming  evidence  for  average,  large-scale  sunward  convection  in 
the  earth's  subpolar  ionosphere  and  inner  magnetosphere  [e.g.,  Stern,  1  0  7  7  ^ . 
We  propose  that  the  most  promising  way  out  of  this  quasi -di  lemma  lies  m  i 
type  of  time-dependent  convection  in  which  plasma  is  released  from  previously 
closed  field  lines  in  sporadic  events,  namely  substorms.  We  do  not  yet 
attempt  to  explain  many  of  the  important  observational  characteristics  of 
substorms  on  this  basis,  but  focus  instead  on  the  proposed  basic  function  of 
the  substonn  in  the  magnetospiier i c-convect i on  process. 


The  Inconsistency 


Consider  slow,  steady,  lossless,  sunward  convection  in  the  Hnd  of  long- 
inagnetotai 1  magnetic-field  configuration  shown  in  figure  2 6,  and,  lor 
simplicity,  assume  isotropic  pressure. 


The  pressure  variation  along  a  drift  trajectory  then  follows  the  law 


P  -  V~5/3  (la) 

where  the  flux-tube  volume  V  is  given  by 

V  /ds/B  (lb) 

and  the  integral  is  along  the  full  length  of  a  closed  magnetic-field  line, 
from  the  southern  ionosphere  to  the  northern  ionosphere.  (See,  e.g.,  Append! » 
of  Hare!  et  al .  (1980a).)  In  a  configuration  like  that  shown  in  figure  26,  we 
define  the  effective  length  of  the  field  line  that  crosses  the  equatorial 
plane  at  x  =  xg  by 


L(xe)  =  V(xe)  B£(xe)  (2) 

where  B ^ ( x^ )  is  the  magnetic  field  in  the  tail  lobe  at  x  -  x(1.  Combining 
(la)  and  (2),  we  obtain  a  prescription  for  the  particle  pressure  in  the 
equatorial  plane  at  x  -  xp : 


»  0(«e)/Bt(xe)r5/3  (3) 


Force  balance  in  the  z-direction  in  the  noon-midnight  meridian  plane  requires 
that 


P  (x  ) 
e  e 


B*(x  ) 
Jl 

'  2iT  ■ 


d z  B  (x  ,z) 
x '  e  ‘ 


B  (x  ,z) 
3x  e  ; 


(4) 


where  z  is  the  z-conrdinale  in  the  tail  lobe,  and  B^(x  )  Bx(xp,z?).  In 

the  derivation  of  (4),  we  have  assumed  B  (x,o,z)  --  B..(x,o,o)  9,  and 

y  A 

l’(x  ,zt)  <<  Pe(xe).  For  xe  <  -20  Rp  ,  the  geometry  is  basically  planar,  the 
magnetic-tension  term  is  small,  and  (4)  reduces  to 


pe(xe)  K  Q)  (S) 


Combining  (3)  and  (6),  we  find  that 
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L(xe)  -x  [B.(xe)]-1/b  (6) 

W\ ■  suggest  (but  cannot  vet  prove)  that  it  is.  impossible  to  construe'  u 
magnet ospheric  magnetic-field  .nodes  tha'  satisfies  (b)  in  the  region  uf  the 
magnetota 1 1  plasma  sheet  and  also  resembles  figure  26  (i.e..  has  a  la  rye 
planetary  field,  a  long  tail  and  a  large  region  of  closed  tail -like  field 

lines).  It  the  tail  is  long  and  the  lobe  is  in  pressure  balance  with  the 

solar  wind,  the  lobe  field  must  decline  slowly  with  distance  down  the  tail, 

for  the  earth,  observations  indicate  that  the  lobe-field  declines  as 

x  (  (e.y.,  Bohannon ,  1968,  Mi hal  ov  and  Sonet t ,  1968).  Equation  (6) 

then  implies  that  L ( xG )  varies  extremely  slowly  with  x(, ,  as  something  between 
the  0.98  and  0.16  power  of  | x J  for  the  case  of  the  earth,  do  feel  that  it  is 
probably  impossible  to  construct  a  magnetospheric  magnetic-field  mode!  in 
which  the  effective  bounce  length  varie-  that  slowly  with  xp.  over  the  entire 
range  of  x0  for  which  the  pressure  relationships  (3)  and  (6)  are  approximate!  v 
valid. 

We  have  tested  several  quant i tat ive  models  of  the  earth’s  magnetospheric 
magnetic  field  for  consistency  with  (1)  and  (4).  Fiyure  27  displays  V(x  ), 
l>  (x  )  and  L(xe)  for  these  models.  Note  that  the  effective  length  L  is 
roughly  linear  in  xp.  This  is  a  much  steeper  dependence  on  xp  than  the 

^  |  /  k 

B,(x  )  dependence,  which  is  also  displayed  for  the  models.  figure  3 
shows  model  values  uf  the  ratio  P  ,(x(,)/Pt.(x,.) ,  where  P  (x  )  is  computed  : ro"i 
(4),  and  P cl ( x e )  is  the  particle  pressure  in  the  equatorial  plane  resulting 
from  adiabatic  convection  from  xe  -  -60  ", 

Pd(xe)  -  Po(-60)  [V(-oO)/V(a  .  ]l'/;  (7) 

Note  that  if  the  two  pressures  balance  a!  xM  -••>0  S',  .  they  disagree  by  at 
least  an  order  of  magnitude  by  x0  -  .  A  flux  tube  shortens 

coris  iderably  as  it  convects  from  B  -60  9r  to  =  -10  9,  ,  which  would  cause  its 
particle  pressure  to  increase  greatly  if  no  particles  are  lost.  There  is  no 
corresponding  major  increase  in  the  confinement  due  to  magnetic  tension  or  the 
magnetic  pressure  of  the  tail  lobe. 

To  put  the  situation  another  way,  suppose  we  adopt  the  01 son-Pf i tzor 
(19/4)  model,  assume  pressure  balance  with  the  lobe  field  at  x(,  -  -60,  and 


also  assume  adiabatic  compression  bet /een  x  =  -60  and  -10,  we  find  that  the 
ratio  6  of  particle  pressure  to  magnetic  pressure  is  ■=  6b  in  the  equatorial 
plane  at  x0  -  -10.  Such  a  high  8  value  is,  as  far  as  we  know,  never  seen  on 
dipole-type  field  lines,  B  ~1  being  much  more  common  in  the  ring  current  and 
the  inner  (dipolar)  part  of  the  plasma  sheet.  Apparently,  flux  tubes  do  not 
convect,  in  an  approximately  lossless  and  adiabatic  way.  from  x(,  «  -60  to  the 
dipole-like  field  region  at  x  =  -10. 

Table  3  lists  the  most  important  simplifying  assumptions  that  we  have 
made  (mostly  implicitly)  in  the  preceding  arguments,  and  comments  on  their 
probable  validity. 

Our  assumptions  are  not  completely  unassai lable.  However,  we  regard  them 
as  adequately  accurate  for  demonstrating  at  least  an  order-of-magn i tude 
pressure  discrepancy.  Given  these  assumptions,  we  can  prove  the  pressure- 
balance  inconsistency  only  for  existing  magnetic-field  models,  not  in  complete 
generality.  However,  given  the  fact  that  all  the  models  exhibit  the  same 
inconsistency  (same  direction,  -  same  magnitude),  and  given  the  discussion 
following  (6),  we  feel  that  the  pressure-bal ance  inconsistency  does  not  result, 
from  simple  defects  in  the  models,  we  believe  the  pressure-balance  inconsis¬ 
tency  to  be  a  general  theoretical  problem  for  a  planetary  magnetosphere  whose* 
particles  convect  planetward,  on  closed  field  lines,  from  far  out  in  a 
magnetotai 1 . 

We  should  mention  a  piece  of  independent,  though  indirect,  observational 
evidence  against  steady,  lossless  convection  in  tho  earth's  magnetotai 1.  The 
observed  latitudinal  distributions  of  region-2  Birkeland  currents  suggest  that 
flux-tube  content  increases  with  distance  out  into  the  tail  [Harel  ot_  al  , 
1980b],  Following  our  intuition  that  tho  pressure-bal ance-i neons l stone v 
argument  is  not  only  valid  but  represents  a  profound  and  general  problem,  we 
devote  the  next  section  to  an  idea  as  to  how  the  earth's  magnetosphere 
eliminates  tho  inconsistency. 

Resolution  of  the  Inconsistency 

Pressures  must  remain  in  balance  in  a  slow  convection  process,  so  a 
slowly  convecting  magnetosphere  must  “"ind  a  way  to  remove  the  pressure-bal ance 
inconsistency  exemplified  by  Figure?*:.  The  most  obvious  way  to  do  this,  while 


retaining  the  picture  of  flux  tubes  t  onvecting  planetward  through  a  lone, 
magnetotail,  is  for  a  large  fraction  ef  particles  to  be  lost  from,  or  de¬ 
energized  on,  converting  flux  tubes.  There  are  various  ways  in  which  this 
might  conceivably  occur.  However,  we  note  that  the  earth's  plasma-sheet 
configuration  is  known  to  vary  greatly  in  time,  particularly  during  substorms, 
and  that  plasma-sheet  flows  are  very  stochastic.  There  is  no  evidence  for 
steady  sunward  flow  in  the  tail  in  quiet  times  (Akasofu,  private  communi¬ 
cation).  We  now  discuss  how,  by  relaxing  the  assumption  of  time  independence, 
we  can  find  a  plausible  picture  of  how  the  magnetosphere  might  release  plasma 
from  middle-plasma-sheet  flux  tubes,  to  allow  these  flux  tubes  to  move  further 
sunward  and  complete  the  convection  cycle. 

The  Ni shida/Hones  substorm  picture  [e.g.,  Nishida  and  Nagayama ,  1  ‘ * 73. 
Hones .  1979,  see  also  Hill,  1973],  is,  along  with  the  Akasofu  picture  [e.g., 
Lin_  et_al.,  1977]  and  the  UCLA  pictun  [e.g.,  McPherron  et  al . ,  1973],  a 
popular  observat i on-based  view  of  plasma-sheet  behavior  in  a  substorm.  I* 
suggests  that  the  release  of  plasma  from  inner-pl asma-sheet  flux  tubes 
constitutes  an  essentia!  part  of  the  substorm  process.  One  possible  scenario 
of  how  the  magnetosphere  might  release  plasma  is  shown  in  figure  29;  this  is 
the  Ni shida/Hones  picture,  slightly  modified  for  the  present  discussion.  The 
top  diagram  shows  an  initial  configuration  corresponding  to  the  average 
magnetic-field  configuration  (presumably  i nner-pl asma-sheet  flux  tubes  like  A 
nave  less  plasma  content  than  the  outer-pl asma-sheet  tubes  like  B  do). 
Convection  proceeds  continuously  in  the  ionosphere.  However,  as  the  out.or- 
plasma-sheet  flux  tube  (B)  convocts  earthward,  the  external  pressure  applied 
to  it  does  not  increase  rapidly  enough  to  cause  it  to  contract  as  much  as 
would  be  required  to  keep  the  overall  magnetic  field  configuration  constant  ’> 
time.  Thus  the  cross-hatched  flux  tube  in  the  third  diagram  has  larger  volume 
and  is  more  stretched  and  tail -like  than  tube  A  in  the  top  diagram,  which  has 
tne  same  ionospheric  foot,  the  same  amount  of  magnetic  flux,  but  less 
plasma.  We  guess  that  this  growth  phase  will  continue,  with  magnetic  field 
lines  in  the  inner  and  middle  plasma  sheet  becoming  more  and  more  stretched 
and  tail-like  with  time,  until  the  plasma  sheet,  becomes  sufficiently  thin, 
and/or  becomes  sufficiently  small  for  some  breakup  mechanism  to  be 

operative.  One  possible  mechanism  is  the  ion-t ear i ng-mode  instability 
discussed  bv  Schindler  (W/4)  whereby  one  or  more  near-earth  neut  ral -1  i  nes 
suddenly  form  (third  and  fourth  diaera  r).  Most  at  the  plasma  from  the  highly 


distended  field  lines  escapes  downtail,  the  closed  field  lines  become  more 
dipolar,  and  the  inner  edge  of  the  plasma  sheet  is  able  to  convect  further 
earthward.  [The  transition  from  tail-like  to  dipole-like  field  lines  is  a 
long  established  characteristic  of  substorms  (Russell  and  McPherron,  1973) , 
and  the  direct  relationship  between  the  location  of  the  inner  edge  of  the 
plasma  sheet  and  flux-tube  plasma  content  was  investigated  quantitative  y  by 
Jaggi  and  Wolf  (1973).]  After  sufficient  plasma  has  been  released  from;  the 
plasma  sheet,  a  recovery  phase  (not  shown  in  Figure  39)  begins  in  whicn  the 
near-earth  neutral  line  probably  moves  tailward  following  the  escaping  plasma, 
and  the  plasma  sheet  is  repopulated  by  the  acceleration  of  particles  in  the 
current  sheet. 

The  theoretical  arguments  advanced  in  this  paper  suggest  that  the 
magnetosphere  will  exhibit  gross  time  variations  even  if  the  solar  wine'  is 
completely  steady,  providing  that  the  steady  wind  drives  nagnetosph>  >  ; 
convection.  Southward  IMF,  which  increases  the  rate  of  convection,  would  m-~d 
to  cause  tne  magnetospheric  disruptions  tc  be  more  frequent  and/or  -v 
intense.  We  suggest  that  the  magnetospheric  substorm  is  a  theo'-et  ica 1 1  y 
necessary  element  of  the  convection  process,  the  element  that  allows  middle- 
plasma-sheet  flux  tubes  to  lose  plasm,  so  that  they  can  move  further  sunward 
and  complete  the  convection  cycle. 
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Table  1.  Model  assumptions. 

A .  Characteristics  of  the  model 

(1)  Region  of  model  =  inner  magnetosphere  and  ionosphere;  L  T  10, 
invariant  latitute  1  70  . 

(2)  0 lson-°f it zer  analytic  magnetic  field  model  and  substorm  curren 
loop  (not  self-consistent). 

(3)  Time  dependent  ionospheric  conductivity  model,  including 
day-night  asymmetry  and  auroral  enhancement. 

(4)  7  •  j  =0  both  in  magnetosphere  and  ionosphere. 

Self-consistent ly  calculated  current  system,  consisting  of 
horiyontal  ionospheric  currents,  magnetospheric  ring  currents, 
and  Birkeland  currents  connecting  the  two. 

(5 )  -  -  -  7  '/in  ionosphere. 

(6)  Isotropic  distribution  of  particle  pitch  angles  in  the 
magnetosphere. 

(7)  (Discrete  particle  energy  spectrum.  Electron  energy  '<1-4  keV; 
ion  energy  v  500  ev  -  60  keV. 

(3)  Electron  loss  via  strong  pitch-angle  scattering. 

(9)  Gradient,  curvature  and  E  x  3  drifts  included.  °o 1 ari cat i on 
currents  excluded. 

( 10) Time- independent  particle  input  at  the  high-L  boundary  of  the 
model.  Maxwellian  input  distributions  are  assumed  with  Tj 
4.5  keV,  Te  1.5  keV  at  L  '  10. 

Not  Included  in  the  Model 
(1)  r ie  Id-al igned  electric  fields. 


8. 
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(2)  Ionospheric  neutral  winds. 

(3)  Pre-existing  ring  current.  All  particles  assumed  to  originate 
at  L  £8. 

(4)  Polar-cap  phenomena.  Solar  wind  and  polar  cap  phenomena  enter 
model  only  through  boundary  conditions. 

(5)  Equatorial  electrojet.  Equatorial  boundary  condition  is 

j  .,=  0  at  2T  0  latitude, 

north-south 
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Table  2.  Summary  of  Major  Comparisons  Between  Model  and  Data  for  the 
Substorm- type  Event  of  19  September  1976 


Simul taneous 


Data  or  Data 

Level  of 

Observationa 1 

Data 

t  rom  similar 

Agreement/ 

Feature 

Source 

Events 

Interpretation 

Our  Reference 

Electric-field 
equatorward  of 
polar-cap 
boundary 

E  detector 
from  S3-2 
(AFGL/ 

Regis 
College 
group  -- 

Simul taneous 

Fairly  good  agreement  but  many 
detailed  discrepancies.  Good 
agreement  on  degree  of  shield¬ 
ing  of  low-latitude  ionosphere 
and  on  dawn-dusk  asymmetry. 

The  most  dramatic  feature  of 

Hard  et  a  1  . 

0  980b) 

W.  J. 

Burke 
et  al . ) 

the  data  -  an  instance  of  large 
(100  mV/m)  electric  field  in 
the  ionospheric  trough  -  was 
well  modeled. 

Birkeland 
currents- 
qual itative 

TRIAD  mag¬ 
netometer 
(APL  group) 

Similar  Events 

Good  Agreement 

Hare!  et  al 

0 98Gb)  " 

overall 

pattern 

e.g.  lijima 
and  Potemra 
(19781 

Birkeland 

currents 

Magneto¬ 
meter  on 
S3-2 

(AFGL/Regis 
Col  lege 

Simultaneous 

Fair  Agreement.  Sense,  loca¬ 
tion  and  average  magnitude 
are  about  right.  Crucial 
theoretical  prediction  of 
downward  current  in  the 

Harel  et  al 
0981 Obj 

group  -- 
F.  J.  Rich 
et  al  .1 

rapid-subauroral-f low- region 
proved  correct.  Magnitudes 
often  off  by  a  factor  of  2 

possibly  due  to  inaccurate 
conduct! vi ty  model  (?).  Also, 
two  interesting  types  of  dis¬ 
crepancies:  (1)  data  show 
Birkeland  current  effects  of 
pre-existing  ring  current, 
which  was  neglected  for 
simplicity  in  these  initial 
models;  (?)  data  show  large 
Birkeland  current  poleward  of 
the  plasma-sheet  inner  edge; 
discrepancy  may  be  due  to  in¬ 
accurate  assumption  of  plasma- 
sheet  energy  distribution  or 
may  indicate  an  additional 
Birkeland  and  current  genera¬ 
tion  mechanism,  in  the  far 
plasma  sheet,  a  mechanism  not. 
included  in  the  model. 

Ring-current  UCSD  Similar  Events  Ring-current  injection  proceeds  Hard  et  al 

particle  fir  modeled  period  of  four  hours  (19R0b) 

detectors  a 'ter  onset.  The  model  produces 

ATS-5,  a  classic  ion-energy  dispersion 

ATS-6  (e.g.  seen  many  times  by  Mcllwain  and 

DeForest  coll abora tors . 

and  Me  I lwa in, 
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Observational 
_  Feature 
Equatorial 
east-west 
electric 
field  at 
L  4 


Dst 


Low-lati¬ 
tude  ground 
magnetograms- 
a symmetric 
ring  current 


AFGL  Magneto- 
grams 


Aurora  1 -zone 
magnetograms 


Table  2  (Cont'd) 


Simul taneous 

Data  or  Data  Level  of 

Data  from  Similar  Agreement/ 

Source _  Events?  _ Interpretation  __  _ Our  Reference 

Whistlers  Similar  Events  Very  good  agreement  between  5piro  et_al. 

(Carpenter  model  and  average  observed  ( 1 980) 

et  al . ,  sub storm  electric  field,  as 

19797  function  of  local  time.  Model 


provides  explanation,  in  terms 
ol  Hall  conductivity  and  day- 
night  asymmetry  in  conductivity, 
of  the  fact  that  the  strong 
westward  electric  fields  are 
observed  in  the  midnight-to- 
dawn  sector,  not  dusk-to-mid- 
n  ijgh  t . _ 


NSSDC /  Simultaneous  Adequate  agreement  with  regard  Section  IV 

M.  Sugiura  to  magnitude  of  decrease  in 

Dst.  Detailed  timing  not 
understood.  Observed  quick 

_  _ recovery  of  Dst  not  understood 

Classic  Similar  Events  Asymmetry  of  model  Birkeland  Section  IV 

observa-  currents  produces  an  asymmetric 

tion  depression  of  the  H-component 


with  strongest  depression  in 
the  dusk-midnight  sector.  This 
is  a  classic  observation  feature 
of  the  early  main  phase  of  a 
_ _  storm. 

P.  Fougere  Simultaneous  Good  agreement  between  model  and  Section  IV 

data  with  regard  to  vertical 
component  of  AB .  The  northward 
and  eastward  AB's  were  smaller 
both  in  the  data  and  in  the 
model ,  where  they  represented 
the  rather  small  net  result  of 
near  cancellations  of  the  mag¬ 
netic  effects  of  large  currents. 

World  Data  Simultaneous  Adequate  agreement  for  Fort  Section  IV 

Center  Churchill  magnetograms.  Poor 

agreement  for  College  magneto¬ 
grams,  due  to  inaccurate  input 
_ _ _ _  a  s  sump t  ions. 


r  b 


Assumption 

No  particle  loss 
during  steady 
convect ion 


Isotropic  particle 
pressure  in  the 
plasma  sheet 


Particles  are  in 
bounce  equi 1 ibrium 


Table  3 

Validity  of  Assumptions 
Comments  on  Validity 


Charge  exchange  is  estimated  to  be  minor  for 
xg  >>  -10,  while  a  large  fraction  of  the  elec¬ 
trons  could  precipitate  for  xe  >  -15  (Kennel  , 

1969).  However,  we  note  that  the  ions  are 
responsible  for  the  bulk  of  the  pressure,  so 
that  electrons  can  be  neglected  for  the 
purposes  of  this  paper. 

Our  drift  calculations  suggest  that  less  than 
about  46%  of  the  particle  pressure  can  be  lost 
due  to  drift  out  the  sides  of  the  tail  between 
xe  =  -60  and  -15.  Even  when  the  losses  mentioned 
here  are  taken  into  account  there  still  exists  an 
an  order  of  magnitude  pressure  discrepancy  at 
x  ^  1?  —15. 

Observations  of  Stiles  et  al .  (1978)  indicate 
that  the  particle  pressure  is  usually  approxi¬ 
mately  isotropic.  We  have  done  pressure-change 
calculations  for  -60  xe  >•  -10  and  various 
equatorial  pitch  angles  assuming  conservation  of 
the  first  two  adiabatic  invariants.  Only  severe 
pitch-angle  anisotropies,  with  pressure-hearing 
particles  confined  quite  close  fn  Hie  center  of 
the  current  sheet,  could  possibly  remove  tv 
pressure  discrepancy.  Such  extreme  an i sot ropies 
seem  inconsistent  with  observation. 

We  have  done  "part icl es- 1 n-a -box"  calcula* ions 
simulating  particles  trapped  in  earthward 
convecting  flux  tubes.  From  these  calculat ions  , 
we  estimate  that  the  bounce-equi 1 ibrium/adiabat ic- 
compression  approximation  gives  an  accurate  enough 
value  (i.e.,  within  a  factor  of  ?)  for  the  pressure 
in  the  equatorial  plane  if  |xej  0.75  ix|M!, 
where  xm  is  the  x- value  of  the  merging  x-line.  We 
estimate  that  these  conditions  usually  hold  lor 
x0  >  -50  or  -of).  In  the  opposite'  case  to  ours, 
where  no  particles  bounce.  Hill  and  Reiff  ( 1 9.’,0a  ,h ) 
have  found  an  inconsistency  somewhat  analogous  *o 
ours  with  the  idea  of  steady  convective  magnet. it 
merging. 


Figure  1.  Corresponding  regions  in  the  ionosphere  and  the  equatorial 
plane;  mapping  is  done  via  the  pre-substorm  magnetic-field 
model.  The  circle  I  in  the  ionosphere  represents  the  as¬ 
sumed  polar  boundary  and  maps  to  curve  I  in  the  equatorial 
plane.  Our  computer  model  applies  to  the  region  equator- 
ward  and  earthward  of  curve  3. 

Figure  2.  Overall  logic  diagram  of  our  program.  Arrows  indicate  flow 
of  information  in  the  program.  Dashed  lines  indicate  fea¬ 
tures  that  we  plan  to  incorporate  in  the  program,  but  are 
not  included  yet.  The  rectangles  at  the  corners  of  the 
central  pentagon  represent  basic  parameters  computed.  In¬ 
put  models  are  indicated  by  rectangles  with  round  corners; 
input  data  are  indicated  by  curly  brackets.  The  program 
cycles  through  the  entire  main  loop  (including  all  the 
rectangular  boxes)  every  time  step  At  (approximately  every 
30  seconds).  The  basic  equations  that  the  computer  uses  or 
solves  are  described  briefly  by  words  or  symbols  next  to  the 
logic-flow  lines. 

Figure  3.  Characteristics  of  the  substorm-type  event  of  19  September 
1976.  The  top  panel  shows  the  H-magnetogram  from  Fort 
Churchill.  In  the  lower  panel,  boxes  represent  polar-cap 
potential  drops  estimated  between  electric-fields  reversals, 
while  error  bars  give  potential  drops  estimated  for  our 
modeling  boundary,  which  coincides  with  the  equatorward  edge 
of  the  region-1  currents. 

Figure  4.  Qualitative  sketch  of  the  ring-current  and  northern-hemi sphere 
Birkeland-current  part  of  our  model  current  system. 


47 


Our  model  current  system  also  includes  southern-hemisphere 
wires  (not  shown).  The  sketch  is,  of  course,  an  over¬ 
simplification  of  the  model  current  system,  which  involves 
about  2700  wires. 

Figure  5.  Theoretical  breakdown  of  contributions  to  the  northward  (x) 
magnetic-field  perturbation  at  Newport,  Washington  (55.08' 
geomagnetic  latitude,  300.01°  geomagnetic  longitude).  The 
magnetic  local  time  was  about  0125  at  substorm  onset  (1000 
UT).  The  code  used  for  the  various  currents  is  as  follows: 

"RC"  =  magnetospheric  ring  current; 

"Rl"  =  region-1  Birkeland  current; 

"EJ"  =  the  part  of  the  northern-hemisphere  auroral  electrojet 
that  lies  poleward  of  our  main  modeling  boundary,  i.e., 
poleward  of  the  equatorward  boundary  of  the  region-1 
currents . 

"N/S"  =  north-south  ionospheric  current. 

"TOTAL"  =  total  theoretical  prediction. 

"E/W"  =  east-west  ionospheric  current  within  the  region  of  our 
main  simulation,  i.e.,  equatorward  of  the  equatorward 
edge  of  the  region-1  currents. 

"R2"  =  region-2  Birkeland  current.  Contributions  from  southern- 
hemisphere  currents  are  shown  in  the  lower  panel. 

Total  AB  and  ring  current  are  shown  with  contributions 
from  northern-hemisphere  currents. 

Figure  6.  Theoretical  breakdown  of  contributions  to  the  eastward  (y)  mag¬ 
netic  perturbations  at  Newport.  The  format  is  the  same  as 
Figure  5. 

Figure  7.  Theoretical  breakdown  of  contributions  to  the  downward  (z)  mag¬ 
netic  perturbations  at  Newport.  The  format  is  similar  to 
Figure  5. 

Figure  8.  Comparison  of  theory  (solid)  and  data  (dashed)  for  the  Newport, 
Washington  magnetogram.  The  lower,  middle  and  upper  boxes  show 
theory  and  data  for  the  northward,  eastward  and  downward 
components,  respectively.  The  zero  levels  are  essentially 
arbitrary  in  each  case,  so  the  curves  were  constrained  to  agree 
at  0900.  "The  dotted  curve  represents  a  typical  quiet  day." 


Figure  9.  Theoretical  breakdown  of  contributions  to  the  northward  (x) 

magnetic  perturbation  at  Camp  Douglas,  Wisconsin  (54.53'  geo¬ 
magnetic  latitude,  333.18  geomagnetic  longitude).  The  mag¬ 
netic  local  time  was  about  0340  at  substorm  onset  (1000  UT). 

The  format  is  similar  to  Figure  5. 

Figure  10.  Theoretical  breakdown  of  contributions  to  the  eastward  (y) 

magnetic  perturbation  at  Camp  Douglas.  The  format  is  similar 
to  Figure  5. 

Figure  11.  Theoretical  breakdown  of  contributions  to  the  downward  (z) 

magnetic  perturbation  at  Camp  Douglas.  The  format  is  similar 
to  Figure  5. 

Figure  12.  Comparison  of  theory  and  data  for  the  Camp  Douglas  magnetogram. 
The  format  is  similar  to  Figure  8. 

Figure  13.  Theoretical  breakdown  of  contributions  to  the  northward  (x) 
magnetic  perturbation  at  Sudbury,  Massachusetts  (53.70  geo¬ 
magnetic  latitude,  357.05  geomagnetic  longitude).  The 
magnetic  local  time  was  about  0515  at  substorm  onset  (1000 
UT). 

Figure  14.  Theoretical  breakdown  of  contributions  to  the  eastward  (y) 
magnetic  perturbation  at  Sudbury.  The  format  is  similar  to 
Figure  5. 

Figure  15.  Theoretical  breakdown  of  contributions  to  the  downward  (z) 
magnetic  perturbation  at  Sudbury.  The  format  is  similar 
to  Figure  5. 

Figure  16.  Comparison  of  theory  and  data  for  the  Sudbury  magnetogram.  The 
format  is  similar  to  Figure  8. 

Figure  17.  Theoretical  breakdown  at  contributions  to  the  northward  (x) 
magnetic  perturbation  at  the  equator  at  1500  Magnetic  Local 
Time.  The  format  is  otherwise  similar  to  Figure  5. 

Figure  18.  Same  as  Figure  17,  but  for  2100  MLT. 

Figure  19.  Same  as  Figure  17,  but  for  0300  MLT. 

Figure  20.  Same  as  Figure  17,  but  for  0900  MLT. 

Figure  21.  Comparison  of  the  average  total  ABx  from  Figures  17-20,  a 

"Theoretical  Dst,"  with  the  observed  Dst.  The  zones  were 
arranged  so  that  the  theoretical  value  at  0900  UT  aqrees  with 
the  observed  Dst  for  0800-0900  UT . 
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Figure  22. 


Figure  23. 


Figure  23. 


Figure  24. 


Figure  25. 
Figure  26. 
figure  27. 


Theoretical  breakdown  of  contributions  to  the  northward  (x) 
magnetic  perturbations  at  Fort  Churchill  (68.7  geomagnetic 
latitude,  322.7°  geomagnetic  longitude).  The  magnetic  local 
time  was  0255  at  substorm  onset  (1000  UT).  The  format  ’ 
similar  to  Figure  5. 

Theoretical  breakdown  of  contributions  to  the  eastward  (y) 
magnetic  perturbations  at  Fort  Churchill.  The  format  is 
similar  to  Figure  5. 

Theoretical  breakdown  of  contributions  to  the  eastward  (\) 
magnetic  perturbations  at  Fort  Churchill.  The  format  ;s  similar 
to  Figure  5. 

Theoretical  breakdown  of  contributions  to  the  downward  (z) 
magnetic  perturbations  at  Fort  Churchill.  The  format  is 
similar  to  Figure  5. 

Comparison  of  theory  and  data  for  the  Fort  Churchill  magneto¬ 
gram.  The  format  is  similar  to  Figure  8. 

View  of  the  assumed  magnetic-field-line  geometry  in  the  earth's 
magnetotail.  The  y  axis  is  directed  out  of  the  page. 

Various  magnetic-field-model  parameters  for  the  magnetotail, 
plotted  against  x,  the  distance  from  the  earth  along  the  earth- 
sun  line,  in  earth  radii.  All  plots  apply  to  the  noon-midnight 
meridian  plane.  "0P74"  and  "B79"  refer  to  the  01 son-Pf i tzer 
(1974)  and  Beard  (1979)  magnetotail  model,  respectively  to  the 
Beard  (1979)  model,  we  have  added  the  earth's  dipole  field. 

(a)  Flux-tube  volume  V,  in  nr/Wh,  as  a  function  of  the  x- 
coordinate  of  the  equatorial  crossing  point. 

(b)  Tail-lobe  magnetic  field  B  in  nanoteslas  (for  z  -  12  Rp). 

(c)  effective  flux-tube  length  L  (solid  lines)  and  also  B 

(dashed  lines),  as  functions  of  the  x-coordinate  of  the 
equatorial  crossing  point.  For  a  system  in  pressure  balance 
with  steady,  lossless  convection  the  L  and  B  curves 

would  have  to  be  parallel.  The  B,-^^  curves  have  been 
normalized  to  agree  with  the  corresponding  t  at  x  -60. 

(d)  The  pressure-balance  parameter  V(x(t)  B  (x  )  ‘.normalized 

to  unity  at  x  =  -60.  This  parameter  would  be  independent 

of  x  m  a  system  that  was  in  pressure  balance, 
e 


ro 

figure  28.  The  ratio  Pd  (x  )/Pl}(*(>)  of  the  equatorial  parti*  !<•  pr»--.-.ur  *• 
resulting  from  lossless,  adiabatic,  earthward  convection  to 
the  equatorial  particle  pressure  required  for  pressure 
balance.  The  ratio  is  normalized  to  unity  at  x0  =  -60  R^. 
This  ratio  would  be  independent  of  xg  for  a  convecting 
system  that  is  in  pressure  balance. 

Figure  29.  Time-dependent  plasma  ejection  from  plasma-sheet  flux  tubes, 
in  the  Nishida/Hones  picture.  The  plasma-sheet  magnetic- 
field  lines  shown  divide  the  plasma  sheet  into  regions  of 
equal  magnetic  flux  (per  unit  local  time),  so  that,  in  the 
top  diagram,  flux  tubes  A  and  B  have  the  same  amount  of 
magnetic  flux.  The  plasma  in  flux  tube  B  gradually  convects 
earthward;  as  time  progresses,  field  lines  of  the  inner 
plasma  sheet  become  more  tail-like  (smaller  G  ),  until  a  new 
near-earth  neutral-line  forms  (third  diagram),  and  most  of 
the  plasma  from  tube  B  is  ejected  down  the  tail. 
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